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We consider linear perturbation equations for long-wavelength scalar metric perturbations in gen- 
eralised gravity, applicable to non-singular cosmological models including a bounce from collapse to 
expansion in the very early universe. We present the general form for the perturbation equations 
which follows from requiring that the inhomogeneous universe on large scales obeys the same local 
equations as the homogeneous Friedmann-Robertson- Walker background cosmology (the separate 
universes approach). In a pseudo- longitudinal gauge this becomes a homogeneous second-order 
differential equation for adiabatic perturbations, which reduces to the usual equation for the lon- 
gitudinal gauge metric perturbation in general relativity with vanishing anisotropic stress. As an 
application we show that the scale-invariant spectrum of perturbations in the longitudinal gauge 
generated during an ekpyrotic collapse are not transferred to the growing mode of adiabatic density 
perturbations in the expanding phase in a simple bounce model. 

PACS numbers: 04.50.Kd, 98.80.Cq larXiv:0801."T667l 



I. INTRODUCTION 

An inflationary expansion in the very early universe can produce the approximately scale-invariant spectrum of 
adiabatic density perturbations observed in the cosmic microwave background [l|, but recently there has been a 
Q-i' debate on whether this spectrum could also be generated in cosmological models where the hot big bang phase is 
~5 ■ preceded by a collapse phase. Among the proposed models are the pre-big bang 0, ekpyrotic and pressureless 

collapse models [H, Q ■ Although many of the models seek to embed the four-dimensional cosmology within a higher 
dimensional theory, most of the quantitative calculations are done within an effective four-dimensional theory which 
O^l ■ reduces to general relativity (GR) at sufficiently low energies or late times. However, some modifications to the 
dynamical equations are required at high energies, close to the bounce, to avoid the singularities present in standard 
cosmological solutions in GR. This has lead to confusion in the literature, as different modifications appear to give 
] different results for the propagation of perturbations through the bounce. In particular there are conflicting claims 
as to whether the dominant mode in the longitudinal gauge metric potential ^E", which has an almost scale- invariant 
spectrum during an ekpyrotic collapse, could in principle give rise to a scale-invariant spectrum in the growing mode 
of adiabatic density perturbations after the bounce 0, 0, S, S EH ■ 
' One approach which has been taken to analyze the transfer of metric fluctuations through the bounce from the 
, collapsing to the expanding phase is to evolve the perturbations through an instantaneous transition along a space- 
like hypersurface using the analog of the Israel junction conditions 14], which describe the matching of two 
solutions of the Einstein equations along a time- like hypersurface. This approach was discussed in Ref. [7||, which 
suggested that generic matching conditions would mix growing and decaying modes across the singular hypersurface, 
I but in Ref. [l^ it was argued that no such mixing occurs for adiabatic perturbations. Another approach is to study 
the evolution of perturbations through a non-singular bounce. This can be done by adding higher order terms in the 
gravitational action [ll,[H;[I3|' but the resulting field equations are complex and model-dependent. 

It is tempting to seek a simpler prescription. In a recent work, Alexander et al [Tlj assumed that the physics 
responsible for the bounce does not affect the form of the evolution equation for the perturbations, and the only 
effect of high-energy modifications at the bounce is to modify the background Hubble rate. They choose a specific 
cosmological background which leads to exactly soluble equations, and in which one finds that the dominant mode of 
^ during an ekpyrotic collapse couples to the growing mode of density perturbations after the bounce. However, we 
shall show that this approach implicitly introduces a non-adiabatic component in the perturbations proportional to 
\1/ on large scales. 

Instead we derive a generalised equation for 5* describing the evolution of adiabatic perturbations in the long- 
wavelength limit for a generalised Friedmann-Robertson- Walker (FRW) background, which reduces to the standard 
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equation when the background is general relativity. This equation is derived by requiring that there exists a long- 
wavelength limit in which the evolution of the perturbed universe is the same as that of the FRW background. That is, 
we assume that the same physics applies to long- wavelength perturbations as applies to the homogeneous background 
cosmology [3, [l^ [2^ . In order to obtain a closed second-order evolution equation for a single metric perturbation we 
use the gauge freedom to work in a pseudo-longitudinal gauge (defined in an Apjjendix) without using the Einstein 
field equations. Adopting the same ansatz for the background Hubble rate as in [ll| we find that the dominant mode 
of \? during an ekpyrotic collapse remains decoupled from the growing mode of density perturbations after the bounce. 

Finally, we will comment on the relation of our work to previous discussions of generalised perturbation equations 
on large scales [20l . [2l| . In particular we identify a generalization of the curvature perturbation on uniform density 
hypersurfaces, ^, and show that this quantity is constant in the long- wavelength limit for adiabatic perturbations, 
even if the gravitational field equations are modified and local conservation of energy is violated. 



II. GENERALISED EQUATIONS FOR THE PERTURBATIONS 



We consider a gravitational theory where homogeneous and isotropic spacetimes obey a Friedmann-type constraint 
equation, determining the expansion rate of comoving worldlines, 9, and an equation for its evolution with respect to 
the proper time, r, along these worldlines, 

6^ = 3f, (1) 

For example, in loop quantum cosmolog y a modified effective Friedmann equation (1) can be derived where / = f{p) 
that leads to a cosmological bounce (e.g. |2^). and in Cardassian models [1^ f{p) cx p+Cp^ have been investigated. In 
both these examples local energy conservation along comoving worldlines then fixes the form of g{p,p). For simplicity 
we neglect spatial curvature, which we expect to be negligible in the early universe in any case. 

In general relativity we then have / = SnGp and g — 8ttG{p + p), where G is Newton's constant, p is the energy 
density and p the isotropic pressure of all types of energy-momentum in the universe. More generally, one can always 
define an effective energy-momentum tensor such that the Einstein tensor G^jy — SnGT^^. From Eqs. ([T]) and ^ we 
can identify an effective density and pressure: 

^ - SttG ' ^ - SttG ■ 

Conservation of the Einstein tensor, V^G^^ = 0, then requires conservation of the effective energy-momentum tensor, 
which implies 

±p<^^^^e(p^« + f^), (4) 



or equivalently, from Eqs. ([T]) and 



^f = -(^9- (5) 



However, in the following we will allow / and g to be arbitrary functions of energy, pressure or other variables. In 
particular we will not assume energy conservation. 

We will assume that the spacetime that evolves into our observable Universe can be described by small perturbations 
about a spatially homogeneous and isotropic spacetime, as in the standard hot big bang model. In particular we 
will consider scalar perturbations about a spatially flat Friedmann-Robertson- Walker (FRW) metric. The linearly 
perturbed line-element is given by [l^, H^l 

ds^ = a^[-(l + 2(t))dr]^ + 2B.,, dx^dr] + {(1 - 2^-)^^ + 2E,,j }dx^dx^] , (6) 

where a{rf) is the background scale factor and a comma denotes a partial derivative with respect to the 3D spatial 
coordinates. In general one should also consider vector and tensor metric perturbations, but at linear order these 
decouple from the scalar modes which describe the primordial density perturbations [25]. 
There is a unit time-like vector field orthogonal to constant-?] spatial hypersurfaces |26i . [27| , 

^^ = -(1-0,-5/), (7) 
a 



3 



whose expansion rate is given by 

= Nf^;^ , (8) 

where a semicolon denotes a covariant 4D derivative. In terms of the scalar metric perturbations we have 

(? = 34r(l-</')"-^' + -V2a, (9) 
a a 

where a prime denotes a derivative with respect to the conformal time rj, denotes the 3D spatial Laplacian and 
V^(T is the shear, with 

a = E'-B. (10) 

At zeroth-order the shear vanishes and the background expansion rate is 9q = 37i/a, where Ti. = a' / a is the conformal 
Hubble parameter. 

For the zeroth-order homogeneous (FRW) background the equations ([T]) and ([2]) can be written as 

= y/o, (11) 

n'-re ^ -yffo. (12) 

Adopting the separate universes viewpoint [l^], the same equations ([T]) and ^ can be used to describe the local 
evolution in an inhomogeneous universe if we work in the long-wavelength limit where physical length scales are 
much larger than the time scale set by the expansion rate, A ^ . This is also sometimes called the ultra-local 
approximation [2^. Of course the expansion time diverges at any cosmological bounce {6 = 0) in the very early 
universe, but the bounce itself must have a characteristic time scale associated with very high energies and thus very 
short time scales. 

We then can apply Eqs. H]) and ^ where we take / — foiv) + ^fiVi^) ^-nd g — go{v) + ^giVi^) and the local 
expansion rate is given, to first-order, by Eq. Q. Neglecting all spatial gradients, we can then write the first-order 
equations in terms of the lapse function 0, its derivative, the curvature perturbation ip and its first and second 
derivatives, 

SHi^' + Hcj)) ^ -yJ/, (13) 

^" -W + nq^' + 2{n' -H^)<j) = Y^g. (14) 

Note that these equations are independent of two of the scalar metric perturbations, B and E in Eq. ([6|), which 
determine the anisotropic shear (jlOp . which vanishes in this long- wavelength limit. 

For adiabatic perturbations on large scales different patches of the inhomogeneous universe follow the same trajec- 
tory in phase space, and the adiabatic perturbations correspond to a perturbation back or forward with respect to the 
homogeneous background along this trajectory (19l] . In this case the hypersurfaces of uniform-^ and uniform- (d^/dr) 
coincide. To first-order this requires Sg/g^ = Sf /Jq. 

More generally, we can write any perturbation 5g as a sum of its adiabatic and non-adiabatic parts, 

Sg = $Sf + Sgn^d , (15) 

JO 

where Jijnad is automatically gauge-invariant. Indeed, if we identify / with an effective density and g — f with an 
effective pressure, then Sgnad — 8ttG[Sp^^ — (pg^ /pg^ )^P^^] — ^^GSp^^. If we assume / = f{p) in Eq. H]) and 
impose energy conservation, so that dp/ dr -\- 9{p + p) = along comoving worldlines, then we would require from 
Eq. ^ that g = {df/dp){p + p) and then Sgnad = {df /dp)Spnad- 

Using Eqs. ^E), and US]), we have from Eqs. ((13]) and (HI]) that 



SHH' -n" -n^ , , nw -h^ 2n'^ - nn" 
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a" 



„ . . I , ^ — (5flnad ■ (16) 

Equation (|16p includes the two gauge-dependent metric perturbations ■0 and (/). However, we still have a choice of 
temporal gauge. If we choose a pseudo-longitudinal gauge in which ijj = (jj = (see Appendix A) we find the following 
evolution equation: 

^ „ mn' - H" - 2^3 , 2H'2 - HH" , , 
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For adiabatic perturbations the right-hand-side vanishes and we have a homogeneous second-order evolution equation 
for '5. We can solve this equation by quadratures to find the general solution [20] 



where C and K are constants of integration. This also has a simpler form [30| : 

(19) 



a'' 



-1 + ^ I a^drj 



Although the differential equation (fT7|) has a singular point when H' — = 0, the solution p9|) is clearly regular 
through a bounce. We show in Appendix A that 5* = CH/a^ is a solution on all scales in the pseudo- longitudinal 
gauge, but the second term is only a solution in the long-wavelength limit. 

A. General relativity limit 

For general relativity we know the full equations for linear perturbations about a FRW metric. The background 
expansion rate obeys the Friedmann constraint and evolution equations, pip and p2p . with fo = SnGpp and go — 
8itG{po +Po), while we have the following perturbation equations, including the spatial gradients, (2^. [27|: 

3H(V'' + rC(j)) - V^ijj - HVV -AnGa'^Sp , (20) 
ip" + 2n^' + n(l)' + i2H' + n^)(l) ^ inGa^iSp+^V^U) , (21) 

o 

where — (l/3)i5ij V^II is the anisotropic stress. It is straightforward to see that Eqs. ifTS]) and (|T4|) reduce to Eqs. ((20|) 
and l[2T|). if we neglect the spatial gradient terms, in general relativity, where 5f = SnGSp and 6g = 8ttG{Sp + 6p). 

General relativity also supplies a momentum constraint equation, and an evolution equation for the anisotropic 
shear potential, a. In the longitudinal gauge, where (7 = 0, this becomes a constraint equation which requires ip = (j) 
in the absence of anisotropic stress, 11 = 0, [2J]. From Eqs. (j20p . pT|) and ([23]) we obtain 



'f" + 3(1 + cDH'f' + [2H' + (1 + 3c2)h2 „ c^V^l^- = AnGa^Spn^^d , (22) 
where here ^' is the curvature perturbation in the longitudinal gauge, and the non-adiabatic pressure perturbation is 

Sp = clSp + Spnad , (23) 

where = p^/ p^ is the adiabatic sound speed of the matter. 

Using Eqs. (fTT|) and (fT2|) it is easy to show that in general relativity Eq. (fT7|) reduces to Eq. ([22|) . neglecting the 
spatial gradients, with Sgnad = SnGSpnud- 

More generally, if we use equations ^ to identify an effective density and pressure on large scales, then one can 
show that our generalised perturbation equation (|17p can be written in a "general relativistic" form 

+ 3(1 + c2°ff)H*' + [2n' + (1 + 3cf^)n^]'9 = AttGo^Sp^JI^ , (24) 

where the effective adiabatic sound speed is 



pf _ nw + re - n " 
pf ~ m{n' - 7e) 



III. APPLICATION TO A SIMPLE BOUNCING COSMOLOGY 

We will now consider the evolution of perturbations in the simple bounce cosmology proposed by Alexander et al in 
Ref. They assumed a barotropic fluid in the universe, with pressure p = wp, and adiabatic density perturbations. 
To model the bounce they assumed a specific ansatz for the background evolution, with the conformal Hubble rate 
given by 

^ = (26) 
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where q = 2/(1 + 3w), for which the scale factor has the solution 

« = «^f4^4)'. (27) 
\vf + VoJ 

where rji and at denote the initial values of the conformal time and scale factor, respectively. At early or late times, 
\r]\ ^ 77o: we recover the general relativistic background evolution 

H - ^ , (28) 
V 

but for \ri\ rjQ we require some unspecified modification of the gravitational field equations to produce the modified 
expansion given in Eq. ((26l) . 

In the general relativistic limit (at early and late times, \ri\ ^ ?7o) we can use either the GR perturbation equa- 
tion or its generalisation, pT)) . to obtain 

^ „ 6(1 + w) , , ^ 

^" + -A ^ 29 

77(1 + 3w) ^ ' 

for adiabatic perturbations in the long- wavelength limit. This has the analytic solution 

^ ^ D± + S±{±7j)-^'' , (30) 

where v = {5 + 3w)/[2(l -I- 3w)], D± and S± are constants of integration and the + and — signs correspond to the 
expanding and contracting phases, respectively. 
In the ekpyrotic limit, w — > 00, we obtain 

■9 = D± + S±{±T^)-K (31) 

Normalising to initial vacuum fluctuations in the ekpyrotic model gives a scale-invariant spectrum of fluctuations in the 
growing mode approaching the bounce, 5" ~ S-{—r])~^ , while the constant mode, Z3_, has a steep blue spectrum 
However, observations constrain the spectrum of perturbations in the adiabatic growing mode in the expanding phase, 
D+. 

If we assume that the perturbations are adiabatic in the long- wavelength limit, that is, Sgnad = in Eq. (|17|) . we 
have 

^' H ^ * H ^ — * = 0. (32) 

In the w 00 limit, substituting in Eq. for Hirj), Eq. ([5^ reduces to 

vI/" + ^^4^^*' = 0. (33) 



^ = Ci + C2 , 2 : (34) 



This has the solution 



where Ci and C2 are constants. Matching the solution in Eq. (I34p to the solution in Eq. (I3ip at both early and late 
times, r] ±00, we get 

D+ = D_ and 5*+ = -S^ , (35) 

therefore we see that there is no mixing between the different modes in the early and late time limits. In particular 
the scale-invariant spectrum during an ekpyrotic collapse, couples only to the decaying mode 5*+ in the expanding 
phase, and the adiabatic growing mode in the expanding phase, D+, has a steep blue spectrum. 

This contrasts with the results presented in (llij . where it was assumed that, although the background evolution 
is modified, Eq. (^5]) . the perturbations are governed by the general relativistic equation (P^ . Neglecting spatial 
gradients and setting Spm^d — they obtain the following solution when w 00 

= B1+B2 arctan f — ) , (36) 
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where Bi and B2 are constants. Matching this solution to the solution in Eq. ([3T|) . at both early and late times, 
r] — > ±00, they obtain 

D+ = D_ + —S- and S+ = -S- . (37) 
Vo 

In this case the mode Z?+ receives a contribution from both modes in the collapsing phase. 

Comparing our Eq. (|17p with Eq. (j22p we see that, by modifying the background evolution but using the unmodified 
GR perturbation equation, the authors of Ref. [ll'| implicitly introduced a non-adiabatic perturbation in the evolution 

Eq. ini), 

^ffnad = K 

Note that Alexander et al [W] show that gradient terms are not important on large scales close to the bounce 
{krjo <^ 1). However, by using the general rclativistic equation they have implicitly introduced a non-adiabatic 
perturbation, (5(7nad, proportional to the longitudinal gauge metric potential and this is why the mode in the 
expanding phase receives a contribution from S- which would not otherwise contribute to the growing mode in the 
expanding phase. 

Mode mixing is certainly possible, indeed it happens in general relativity, but only on finite scales due to V^^* ^ 
(see e.g., [3l|). But mixing occurs due to spatial gradients and hence is suppressed on large scales. Of course we 
cannot say that there is no theory in which a non-adiabatic perturbation of the form given in Eq. psp could occur, 
but one would need a physical model for the appearance of non-adiabatic source terms unsuppressed on large scales 
(e.g., a non-adiabatic perturbation of a multi-component system). We have shown that by assuming only adiabatic 
perturbations on large scales, derived from linear perturbations about the background solutions to Eqs. ([T]) and ([2]), 
no such mixing occurs. 



37^(1 



W -TP 



(38) 



IV. CONCLUSIONS 



In this paper we have derived a generalised equation for the evolution of scalar metric perturbations in the long- 
wavelength limit, for a generalised Friedmann-Robertson- Walker (FRW) background, which reduces to the standard 
equation for the longitudinal metric perturbation i n g eneral relativity. 

Our results are consistent with previous work ^lol . |2]| which pointed out that the curvature perturbation on uniform- 
density hypersurfaces [l^, , 

C = -^-^5p, (39) 

is conserved for adiabatic perturbations on large scales assuming only local conservation of energy (see also [3^. [ssl. [3i| ) . 
We can define a generalization 

C/ ^ - '^Sf (40) 
/o 

which is the gauge-invariant definition of the curvature perturbation, —tp, on uniform-expansion hypersurfaces, where 
5f — 0. Using Eqs. (fTTj) and we can write 

C/ = ^^7-^ [wv-' + {n^ - n')^j + n^<p\ . (41) 

In general relativity the uniform-density and uniform-expansion hypersurfaces coincide in the long-wavelength limit 
and hence Cf = C- Using Eqs. ([TH)) and for the evolution of perturbations on large scales we obtain 

C/ = j^, _ ^2 y^ffnad ■ (42) 

We see that C/ is constant in the large-scale limit for modified gravitational field equations, even allowing for non- 
conservation of energy, if the perturbations are adiabatic, i.e., ^^nad = in Eq. (jlSp. 

We see that the growing mode solution (in an expanding universe) for the pseudo-longitudinal gauge perturbation, 
(X K in Eq. (|18p. corresponds to C/ = K, where if is a constant of integration. The decaying mode = CTL/a? 
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does not contribute to the curvature perturbation Q and is the residual gauge mode in the definition of the pseudo- 
longitudinal gauge (see Appendix A) . The other physical degree of freedom in the pseudo- longitudinal gauge resides 
in the anisotropic shear potential, a, which does not affect Q or the large-scale evolution equation (fT7| in the long- 
wavelength limit. 

The same solution (jlSp was given by Bertschinger [20| in modified gravity for the longitudinal gauge metric per- 
turbation assuming that the trace-free spatial part of the Einstein tensor vanishes (equivalent to assuming vanishing 
anisotropic stress in GR). We have shown that it is possible to use the solution ([T5)) for adiabatic cosmological per- 
turbations on large scales in a pseudo-longitudinal gauge without making any assumption about the gravity theory, 
energy conservation or anisotropic stress. The GR solution for the shear in the pseudo-longitudinal gauge in the 
absence of anisotropic stress (e.g., in a conventional hot big bang phase) is given by a — C/a^, where C is a constant 
of integration. In this case we can perform a further gauge transformation drji = a [sec Eq. (IA3p ] to obtain the usual 
longitudinal gauge curvature perturbation, "i/g — + CTL/a? , which shows that in the GR limit the longitudinal and 
pseudo-longitudinal solutions differ only in the decaying mode, oc Ti/a?. 

We applied our formalism to adiabatic perturbations in a simple cosmological bounce model, assuming a specific 
ansatz for the background evolution introduced in Ref. [lH. We took the limit where the equation of state of the 
matter present in the universe w = p/p goes to infinity, corresponding to an ekyprotic model. We conclude that the 
dominant mode of the curvature perturbation after the bounce does not receive a contribution from the growing mode 
in the collapse phase. This result is in contrast with the result presented in Ref. [Till, where they have assumed that, 
although the background evolution is modified, the perturbations are governed by the same equation as in general 
relativity. We have shown that such an approach is equivalent to introducing a non-adiabatic component proportional 
to on large scales. 

Since it is not possible to unambiguously separate perturbations from the background in the large scale limit 
where spatial gradients are negligible, we argue that if the background equations are modified then the perturbation 
equations must also be modified in the large scale limit. And Eq. p?]) . not Eq. (|22p . is the consistent generalisation 
of the evolution equation for long- wavelength linear perturbations given a modified background solution for 'H{r]). 

The general nature of the perturbation equation (|17p and solution (|19p. derived only assuming generalised gravity 
field equations of the form given in ([1]) and implies that the growing mode of metric perturbations (j3ip in an 
ekpyrotic phase will only contribute to the decaying mode in an expanding phase, independently of the details of the 
bounce or the expanding phase. 
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APPENDIX A: PSEUDO-LONGITUDINAL GAUGE 



The metric perturbations and ip given in Eq. ^ and a defined in Eq. pH)) are gauge-dependent. Under a change 
of time coordinate -q ^ rj + Srj we have [13, [13] 

(j, 4) = (l)-ndr]-5r]' , (Al) 

ip = ^j + n6Tj, (A2) 

a a — a — Sr/ . (A3) 



Therefore we can always choose 



Sr]' + 2n5r] = (f> - i/j (A4) 



such that from Eqs. (|Aip and (|A2p we have ip — 4> = '9 . This fixes the gauge up to an arbitrary constant of integration 
C(x), which is similar to what occurs in the synchronous gauge. Perturbations in this pseudo-longitudinal gauge are 
given in terms of metric perturbations in an arbitrary gauge by substituting Eq. (|A4p in either Eq. (|Aip or (|A2p . to 
obtain 

^^, = ^ + li ( a\(j) -i})df]. (A5) 
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We see that we can write the general sohition for vj/ as ^+{ri) + 'i'-{ri), where ^'-(r/) = CT-i/a? decays in an expanding 
universe for w > —1. 

In general relativity we have a field equation coming from the trace-free part of the Einstein tensor which requires 

^ - = STrGa^n -a' - 2ha , (A6) 

where the scalar part of the anisotropic stress is given by ViVjll — {l/3)6ijV^Il. In the longitudinal or Newtonian 
gauge [131 we have = by definition, and hence ijje — 4>e. in the absence of anisotropic stress. More generally, we 
can define a pseudo-longitudinal gauge in which ip = (j) and the shear potential is given by 

ct' + 2ha ^ 8nGa^U . (A7) 

Therefore, in the absence of anisotropic stress, 11 = 0, we have a oc a^^. 

In generalised gravity we may not have a field equation of the same form as Eq. (jA6p , so we cannot in general identify 
the pseudo-longitudinal gauge ip = (p with a zero-shear gauge. If we assume that there exists a general relativistic 
limit in which the pseudo-longitudinal gauge reduces to the usual longitudinal gauge (in the absence of anisotropic 
stress) then this fixes the otherwise arbitrary constant of integration, C(x), and fixes the gauge throughout. 
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